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Abstract. The channeling of the ion recoiling after a collision with a WIMP changes the ionization 
signal in direct detection experiments, producing a larger signal than otherwise expected. We give 
estimates of the fraction of channeled recoiling ions in Nal (Tl), Si and Ge crystals using analytic 
models produced since the 1960's and 70's to describe channeling and blocking effects. We find that 
the channeling fraction of recoiling lattice nuclei is smaller than that of ions that are injected into 
the crystal and that it is strongly temperature dependent. 
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INTRODUCTION 

Channeling and blocking effects in crystals refer to the orientation dependence of 
charged ion penetration in crystals. In the "channeling effect" ions incident upon a crys- 
tal along symmetry axes and planes suffer a series of small-angle scatterings that main- 
tain them in the open "channels" between the rows or planes of lattice atoms and thus 
penetrate much further into the crystal than in other directions. Channeled incoming 
ions do not get close to lattice sites, where they would be deflected at large angles. The 
"blocking effect" consists in a reduction of the flux of ions originating in lattice sites 
along symmetry axes and planes, creating what is called a "blocking dip" in the flux 
of ions exiting from a thin enough crystal as a function of the angle with respect to a 
symmetry axis or plane. These effects were first observed in the 1960's and are used 
in crystallography, in the study of lattice disorder, ion implantation, surfaces, interfaces 
and epitaxial layers, in measurements of short nuclear lifetimes etc. In particular, avoid- 
ing channeling is essential in the manufacturing of semiconductor devices, since ion 
implantation at a controlled depth is the primary technique. 

Ion channeling in Nal (Tl) was first observed in 1973 by Altman et al. [1]. They ob- 
served that channeled ions produce more scintillation light because they lose most of 
their energy via electronic stopping rather than nuclear stopping. The potential impor- 
tance of the channeling effect for direct dark matter detection was first pointed out by 
Sekiya et al. [2] and subsequently for Nal (Tl) by Drobyshevski [3] and by the DAMA 
collaboration [4]. When Na or I ions recoiling after a collision with a dark matter WIMP 
(Weakly Interacting Massive Particle) move along crystal axes and planes, their quench- 
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ing factor is approximately Q = I instead of (2Na = 0-3 and Qi = 0.09, since they give 
their energy to electrons. 

Ion channeling in crystals could give rise to a daily modulation due to the preferred 
direction of the dark matter flux arriving on the Earth. Earth's daily rotation naturally 
changes the direction of the "WIMP wind" with respect to the crystal axes, which 
produces a daily modulation in the measured recoil energy (equivalent to a modulation 
of the factor Q). Sekiya et al. [2] pointed this out for stilbene crystals, and Avignone, 
Creswick, and Nussinov [5] for Nal (Tl) crystals, although the available estimates of the 
strength of these daily modulations are somewhat simplistic. My collaborators, Graciela 
Gelmini and Paolo Gondolo, and I [6] present here our analytic calculations of the 
channeling fraction in Nal (Tl), Si and Ge crystals. 



Our calculation is based on the classical analytic models developed in the 1960's and 
70's, in particular by Lindhard [7-14]. We use the continuum string and plane model, 
in which the screened Thomas-Fermi potential is averaged over a direction parallel to a 
row or a plane. This averaged potential U is considered to be uniformly smeared along 
the row or plane of atoms, which is a good approximation if the propagating ion interacts 
with many lattice atoms in the row or plane by a correlated series of many consecutive 
glancing collisions with lattice atoms. Just one row or plane is considered in this model. 
Lindhard proved that for an ion propagating with kinetic energy E, and for small angle 
between the ion's trajectory and the atomic row (or plane) in the direction perpendicular 
to the row (or plane), the so called "transverse energy", Ej_ = E sin^ + ~ £0^ + U 
is conserved. 

The conservation of the transverse energy provides a definition of the minimum 

distance of approach to the string (or plane) of atoms, Pmin^ at which the trajectory of 
the ion makes a zero angle with the string (or plane), and also of the angle i// at which 
the ion exits from the string (or plane), i.e. far away from it where U ~ 0. In reality the 
furthest position from a string or plane of atoms is the middle of the channel, whose 
width we call Jch- Thus, 



Channeling requires that p^in > pdE), where pciE) is the smallest possible minimum 
distance of approach of the propagating ion with the row (or plane) for a given energy 
E. Since the potential U{p) decreases monotonically with increasing p, ?7(Pmin) < 
U{pc{E)). Using Eq. 1, this can be further translated into an upper bound on Ex_ and 
on ^f, the angle the ion makes with the string far away from it. 



MODELING OF CHANNELING 



E^ = U (pmin) =Ey^ + U ( Jch/2) . 



(1) 




(2) 



y/ciE) is the maximum angle the ion can make with the string far away from it (i.e. 

in the middle of the channel) if the ion is channeled. At low enough E, pc{E) becomes 
close to dch/2, and the critical angle yfc{E) goes to zero. This means that there is a 
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FIGURE 1. Comparison of theoretical (black lines) temperature corrected critical angles (with c = 2) 
and measured critical angles extracted from thermal wave measurements [14] (green, or gray if color not 
available, dots, joined by straight lines to guide the eye) as a function of the energy of (a) B ions and (b) 
P ions propagating in a Si crystal at T=20 °C, for the indicated axial and planar channels. 



minimum energy below which channeling cannot happen, even for ions moving initially 
in the middle of the channel. 

So far we have been considering static strings and planes, but the atoms in a crystal 
are actually vibrating. We use the Debye model, and take into account thermal effects 
in the crystal through a modification of the critical distances which was found originally 
by Morgan and Van Vliet [9] and later by Hobler [14] to provide good agreement with 
simulations and data. It consists of taking the temperature corrected critical distance 
Pc{T) to be, 

p,(r) = ^pm + [cuiiT)]^. (3) 

where ui is the one dimensional rms vibration amplitude of the atoms in a crystal, and 
the factor c in different references is a number between 1 and 2 [9, 10]. As shown in 
Fig. 1, with this formalism and using c = 2 we fit relatively well the critical angles 
measured at room temperature for B and P ions in a Si crystal in several channels, for 
energies between 20 keV and 600 keV that Hobler [14] extracted from thermal wave 
measurements. 

As an example, the static (c = 0) axial and planar critical distances are presented 
in Fig. 2(a) for the 100 channel in Nal crystal, together with the amplitude of thermal 
vibration ui at 20 °C, and the Thomas-Fermi screening distances for Na and I ions. 
Fig. 2(b) shows the temperature corrected axial and planar critical angles at 20 °C (with 
c = 1) for the same channel as functions of energy of the traveling Na and I ions. 

CHANNELING OF INCIDENT PARTICLES 

The channeling of ions in a crystal depends not only on the angle their initial trajectory 
makes with strings or planes in the crystal, but also on their initial position. Ions which 
start their motion close to the center of a channel, far from a string or plane, where they 
make an angle i//, are channeled if the angle is smaller than a critical angle (as explain 
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FIGURE 2. (a) Static critical distances of approach and ui at 20 °C and (b) critical channeling angles 
at 20 °C with c = 1 as a function of the energy of propagating Na (green/gray) and I (black) ions in the 
<100> axial and { 100} planar channels. The screening radii shown as vertical lines are ONa — 0.00878 
nm and o/ = 0.01 15 nm [6]. 



earlier) and are not channeled otherwise. Particles which start their motion in the middle 
of a channel (as opposed to a lattice site) must be incident upon the crystal. 

No data or simulations of Na and I ions propagating in a NaT crystal is available at low 
energies. We show that to a good approximation we can use analytic calculations and 
reproduce the channeling fraction in NaT presented in Ref . [4] . Temperature corrections 
are neglected and we use a static lattice (similar to the approach of Ref. [4]). For an 
incident angle i// with respect to each of the channels and an ion energy E, the fraction 
;^inc(£', y^) of channeled incident ions for axial and planar channels is Xinc = 1 if is 
smaller than the critical angle for the corresponding channel and zero otherwise. 

To find the total fraction of channeled incident nuclei, we average Xinc over the in- 
cident direction. The integrals cannot be solved analytically, so we integrated numeri- 
cally by performing a Riemann sum once the sphere of directions has been divided us- 
ing the Hierarchical Equal Area iso-Latitude Pixelization (HEALPix) method [15]. The 
HEALPix method uses an algorithm to deal with the pixelization of data on a sphere, and 
it is useful to compute integrals over direction by dividing the surface of a sphere into 
many pixels, computing the integrand at each pixel (i.e. each direction, see Fig. 3(a)), 
and finally summing up the values of all pixels over the sphere. 

A channeled ion can be pushed out of a channel by an interaction with an impurity 
such as the atoms of Tl in Nal (Tl). Here we will simply assume that if a channeled 
ion interacts with a Tl atom it becomes dechanneled and thus it does not contribute to 
the fully channeled fraction any longer. We thus neglect the possibility that after the 
interaction the ion may reenter into a channel, either the same or another. Fig. 3(a) 
shows the axial and planar channels of the NaT crystal for incoming Na ions with an 
energy of 50 keV. We include here only the channels with lower crystallographic indices 
100, 110 and 111. The fraction of channeled incoming Na and I ions after including 
dechanneling in the way just described is shown in Fig. 3(b). Our results agree well with 
those published by the DAM A collaboration [4], also included in the figure. 
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FIGURE 3. (a) Channeling fraction for a 50 keV Na ion in different directions plotted on a sphere 
using the HEALPix pixelization: probability equal to one in red, and probability equal to zero in blue, (b) 
Fraction of channeled incident 1 (black) and Na (green/gray) ions as a function of their incident energy 
E with the static lattice without (dot dashed lines) and with (solid lines) dechanneling due to interactions 
with Tl impurities. The results of DAMA are also included (dashed lines). 



The recoiling nuclei start initially from lattice sites (or very close to them), thus blocking 
effects are important. In fact, as argued originally by Lindhard [7], in a perfect lattice 
and in the absence of energy-loss processes the probability that a particle starting from 
a lattice site is channeled would be zero. The argument uses statistical mechanics in 
which the probability of particle paths related by time-reversal is the same. Thus the 
probability of an incoming ion to have a particular path within the crystal is the same 
as the probability of the same ion to move backwards along the same path [16]. This is 
what Lindhard called the "Rule of Reversibility". Using this rule, since the probability of 
an incoming channeled ion to get very close to a lattice site is zero, the probability of the 
same ion to move in the time-reversed path, starting at a nuclear site and ending inside a 
channel, is zero too. However, any departure of the actual lattice from a perfect lattice, 
for example due to vibrations of the atoms in the lattice, would violate the conditions 
of this argument and allow for some of the recoiling lattice nuclei to be channeled, as 
already understood in the 70's [12, 17]. We now estimate the effect using the formalism 
presented so far. 

For axial channels, the probability distribution function g{r) of the perpendicular 
distance to the row of the colliding atom due to thermal vibrations can be represented 
by a two-dimensional Gaussian. The channeled fraction of nuclei with recoil energy E 
making an initial angle with respect to the axis is given by the fraction of nuclei which 
can be found at a distance r larger than minimum distance rj^niin from the row at the 



CHANNELING OF RECOILING LATTICE NUCLEI 
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FIGURE 4. Channeling fractions at T=293 K for Na (solid lines) and I (dashed lines) ions for c — I 
(black) and c = 2 (green/gray) cases (a) without and (b) with dechanneling included. 
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Figure 4 shows what we consider to be our main predictions for the range expected 
as an upper limit to the channeling fraction in Nal (Tl) at 239 K for two different 
assumptions for the effect of thermal vibrations in the lattice, which depend on the 
value of the parameter c used in the temperature corrected critical distances of approach. 
Dechanneling is ignored in Fig. 4(a) and taken into account in Fig. 4(b). 

Our results for the geometric total channeling fraction for Si ions propagating in a 
Si crystal and Ge ions propagating in a Ge crystal at different temperatures are shown 
in Figs. 5 and 6 for the two cases of c = 1 and c = 2, respectively. Please note that we 
have not considered the possibility of dechanneling of initially channeled ions, due to 
imperfections in Si and Ge crystals. Any mechanism of dechanneling will decrease the 
fractions obtained here for Si and Ge. 



We have studied the channeling of ions recoiling after collisions with WIMPs within 
Nal (Tl), Si and Ge crystals. Channeled ions move within the crystal along symmetry 
axes and planes and suffer a series of small-angle scatterings that maintain them in the 
open "channels" between the rows or planes of lattice atoms and thus penetrate much 
further into the crystal than in other directions. Ions which start their motion close to 
the center of a channel, at an initial angle i//, are channeled if the initial angle is smaller 
than the critical angle in Eq. 2, and are not channeled otherwise. We have found that the 
channeling of lattice ions recoiling after a collision with a WIMP is very different from 
the case of incident ions, and that the channeling fraction is smaller. 

As argued originally by Lindhard [7], in a perfect lattice and in the absence of energy- 
loss processes, the probability that a particle starting from a lattice site is channeled 
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CONCLUSIONS 
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FIGURE 5. Channeling fractions of (a) Si and (b) Ge recoils in a Si and a Ge crystal respectively, as a 
function of the ion energy for temperatures T=900 °C (orange or medium gray), T=600 °C (green or light 
gray), 293 K (black), and 44 mK (blue or dark gray) in the approximation of c = 1 . 
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FIGURE 6. Same as Fig. 5 but fore = 2. 



would be zero. However, due to vibrations in the crystal, the atom that interacts with a 
WIMP may be displaced from its position in a perfect lattice, and there is a non-zero 
probability of channeling (as given in Eq. 4). 

As we see in Fig. 4(a) for Nal, without including dechanneling, the channeling frac- 
tion is never larger than 5% and the maximum happens at lOO's of keV. This maximum 
occurs because the critical distances decrease with the ion energy E, making channel- 
ing more probable, and the critical angles also decrease with E, making channeling less 
probable. The simple extreme model of dechanneling we used for Nal predicts much 
smaller fractions, at most in the 0.1% level, with the maximum shifted to small ener- 
gies, less than 10 keV (see Fig. 4(b)). This reduction may eventually prove to be too 
extreme and at present we do not have a better formalism to model dechanneling. With 
the simple model of dechanneling we used for Nal, we could reproduce the channeling 
fractions computed by the DAMA collaboration (which, however, apply to ions which 
start their motion close to the middle of a channel and not to the case of WIMP direct 



detection). 

If the values found by Hobler [14], and also by us above (see Fig. 1), to reproduce 
measured channeling angles in B and P propagating in Si apply also to the propagation of 
Si ions in Si, then the case of c = 2 should be chosen and the channeling fractions would 
never be larger than 0.3%. Moreover, increasing the temperature of a crystal usually 
increases the fraction of channeled recoiling ions (see Fig. 5), but not always. Sometimes 
the opposite happens (see Fig. 6). The c = 1 choice leads to channeling fractions close 
to 1% for Si and Ge. 

The analytical approach used here can successfully describe qualitative features of 
the channeling and blocking effects, but should be complemented by data fitting of 
parameters and by simulations to obtain a good quantitative description too. Thus our 
results should in the last instance be checked by using some of the many sophisticated 
Monte Carlo simulation programs. 
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